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Abstract 

The  detection  of  a  marine  target  with  the  aid  of  an  active  sonar 
necessitates  a  probabilistic  model  of  the  phenomena,  and  detection  improves 
according  to  how  well  the  model  takes  account  of  the  various  random  and 
deterministic  physical  phenomena.  The  signal  wave  Z(t)  received  from  a 
point  source  x  of  the  target  is  described  by  a  deterministic  scalar  or 
vector  function  G(t,x).  Random  reflection  is  characterized  by  a  scalar 
multiplying  factor  C(x),  randomly  distributed  according  to  a  mixture  of 
central  Gaussian  laws.  Also  considered  is  the  model  consisting  of  a  linear 
process  of  the  type  Z(t)  =  /G(t,x)dL(x)  where  L  is  a  Levy  process 
(possibly  non-stationary)  which  describes  the  joint  distribution  of  the 
point  sources  and  their  intensities.  The  problem  consists  of  identifying 
the  model  most  closely  compatible  with  the  experimental  results. 


Introduction 


The  better  the  model  of  the  reflected  sonar  wave  from  a  marine  target 
the  more  efficient  is  its  detection.  The  model  must  be  true  and  precise 

It  is  essential  to  understand  the  physics  of  the  phenomena  in  order  to 
specify  the  mathematical  model.  The  parameters  of  this  model  must  be  determined 
from  experimental  results  and  this  procedure  must  lead  to  formulae  exploitable 
by  computer.  More  precisely  the  statistical  law  of  the  observations  is 
obtained  by  means  of  its  characteristic  function. 

This  works  extends  a  more  theoretical  study  "Identification  of  a  class  of 
filtered  Poisson  processes”  by  the  author  and  Antonio  Gualtierroti . 

The  target  consists  of  points  shining  under  the  effect  of  the  incident 
wave  from  the  active  sonar;  certain  points  reflect  sufficiently  in  the 
direction  of  the  receivers. 

Paragraph  I  is  devoted  to  an  elementary  description  of  the  physical 

phenomena.  A  geometrical  (ray)  approximation  is  used  to  calculate  the  phase 

changes  in  the  emitted  wave  (supposed  sinusoidal) .  Thus  we  find  a  response 

(1.1)  P(t)  =  C(x)G(t,x)  received  at  instant  t  from  the  point  x  of  the  target. 

The  number  N  of  point  sources,  their  positions  x  in  the  target,  and  the 

coefficients  C(x)  are  all  random.  A  process  with  distribution  £  C(x.)6  (x) 

i  i 

(a  Levy  process)  may  describe  these  random  phenomena.  We  suppose  that  effects 
due  to  the  different  sources  are  additive.  This  corresponds  to  the  hypothesis 
of  small  movements  in  mechanics. 

Paragraph  II  defines  the  probability  laws  of  N,  X  and  C.  Three  distinct 
situations  are  described  as  are  the  experimental  results  which  determine  the 
most  realistic  and  the  most  tractable  model. 

In  Paragraph  III  the  variance  and  the  characteristic  function  of  the 
observation  are  treated  case  by  case.  Finally  the  identification  of  the 
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probability  laws  is  presented  in  the  last  paragraph.  It  seems  that  the 
function  G(t,x)  may  not  be  statistically  estimated  from  the  observations 
alone;  a  knowledge  of  the  physics  of  the  phenomena  must  furnish  the  function 
G(t,x),  otherwise  some  very  restrictive  hypotheses  on  the  nature  of  G  would 
have  to  be  imposed. 

The  first  step  of  the  identification  is  to  specify  the  deterministic 
function  G(t,x)  using  physical  considerations. 

Secondly  it  must  be  established  from  the  normalized  observations  whether 
the  law  is  symmetric.  In  the  non-symmetric  case  only  the  model  no.  2  may  be 
applied . 

2 

If  symmetric  (the  most  plausible  case)  the  characteristic  function  of  v=u 

is  tested  to  see  if  it  is  the  Laplace  transform  of  a  probability--if  so,  models 

1  and  3  may  be  used  (cf.  IV. 12  or  IV. 14).  If  in  the  more  restricted  case  the 

2 

second  characteristic  function  of  v=u  is  the  Laplace  transform  of  a 
probability  (cf.  IV. 13)  then  the  second  model  may  be  used  with  a  symmetric 
Levy  process  composed  of  spherically  invariants  defined  by  the  probability 
A.y  (see  formula  III. 9). 

Proposition  IV. 3  gives  a  criterion  for  compatibility  of  a  given  model  with 
the  experimental  results  and  to  know  whether  a  moment  method  or  a  rapid 
Fourier  transform  method  may  be  used  to  find  A.y. 
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I .  Description  of  the  phenomena 

The  submarine  (target)  is  represented  by  a  segment  M  =  [-&,+£]  whose 
centre  is  at  a  distance  R  from  the  active  sonar  S.  In  the  plane  defined  by 
the  sonar  and  the  target,  the  disposition  is  as  in  Figure  1. 


Figure  1. 


Active  Sonar 


Target 

00  =  x  e  [-&,£] 


The  distance  R  between  sonar  and  target  is  supposed  to  be  large  with 


respect  to  the  dimension  i  and  if  —  <<  1  the  emitted  accoustic  wave  is  supposed 
plane.  (An  approximation  by  spherical  waves  would  modify  the  deterministic 
function  G(t,x)  but  not  the  probabilistic  modelling  of  the  phenomena). 

With  this  approximation,  let  us  follow  the  emitted  wave  until  it  is 
received  again. 


Emission :  P(t)  =  Pq  cos  (cot  +  <j)). 

At  point  T  of  the  target: 

-y(R  +  x  sin  0) 


where  co  =  angular  velocity 
<j>  =  phase  at  t=0 


P(t)  =  Pr 


cos 


CO  (t 


R  +  x  sin  0 


o  (R  +  x  sin  0) 

Distance  SN  has  been  approximated  by  R  +  x  sin 
Energy  losses  involve  a  factor  e  ^  per  unit  length. 
The  celerity  is  denoted  -c. 


-)  + 
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Reemission  at  point  T: 

Only  a  part  of  the  energy  is  retransmitted  in  the  direction  S.  We  intro 


duce  a  random  multiplicative  factor  C(x):  thus 

-y(R  +  x  sin  0) 
■  PoC«  e(R  +  x  sin  9T 


COS 


R  +  x  sin  0. 
o)(t - - - j  + 


Reception  at  point  S: 

The  wave  emitted  from  point  T  returns  to  points 


-2y(R  +  x  sin  0) 


cos 


.  2[R  +  x  sin  ). 

“Ct - c - } 


P(t)  =  P  C(x)  -7-=- - : - TT-5- 

v  J  o  K  1  (R  +  x  sm  0)2 

For  simplicity  the  receivers  are  at  the  same  point  S  as  the  emitters. 

This  simplified  analysis  introduces  the  function  C(t,x)  defined  by  the 
equations 

(1.1)  P(t)  =  C(x)G(t,x) 


(1.2)  G(t,x)  =  Pr 


-2y(R  +  x  sin  0) 


cos 


cot  -  2cox  sin  0  - 


2uR 


*♦1 


o  (R  +  x  sin  0j^~ 

It  is  possible  to  introduce  a  function  G(l-x)  which  necessitates  new  approxi- 

% 

mations  (cf.  Westcott) .  Again  for  «  1  <=>  x  sin  0  <<  R,  and  putting 
2<oR 


* 


♦  - 


we  have 


G(t,x) 


, 

0  R2 


cos 


2a)X  sin  0 

U)t - 3 - + 


The  position  x  of  T  may  be  replaced  by  a  temporal  equivalent  by  putting 
2x  sin'  0 


t(x)  = 


from  which 


-2yR 

G(t,n)  =  G (t-t (x) )  ^P 

R 


cos 


co(t-t(x))  +  4» 


The  sources  instead  of  being  situated  at  points  x  of  the  target  would  be  at 
dates  t  =  t(x).  These  last  approximations  do  not  in  fact  seem  to  be 
particularly  useful  in  what  follows. 
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II .  Probabilistic  model 

Certain  points  T\  at  positions  x^  of  the  target  emit  in  direction  T\S. 

We  suppose  the  effects  of  these  different  points  to  be  additive  and  we  take 
the  model 

(11. 1)  Z(t)  =  l  C(x,)G(t,x  )  =  /  G(t,x)dL(x) 

3  J  J 

where  the  process  L  is  defined  on  the  target  (that  is  for  x  e  M  =  [-£,+£]). 

(11. 2)  L (x)  =  l  C(x  )  6  (x) . 

j  3  j 

Stationarity  of  L  has  not  been  supposed  a  priori. 

Case  no.  1. 

Each  observation  at  S  corresponds  to  one  bright  point  T  at  x  on  the  target. 
The  distribution  of  x  is  given  by  a  probability  v(dx).  In  this  case 

(11. 3)  Z(t)  =  C(x)G(t,x)  . 

It  remains  to  specify  the  law  of  C  when  x  is  given.  So  as  not  to  be 
limited  to  a  simple  centred  Gaussian  law,  we  may  suppose  C  to  be  spherically 
invariant  with  associated  probability  (x,da) .  Thus 

(11. 4)  E(exp  iu  C(x))  =  J^+e  a  ^2y(x,da)  . 

More  precisely  is  a  transition  probability.  In  this  way  the  law  of  C  is 
centred  and  symmetric  but  the  distribution  y(x,*)  of  the  standard  deviations 
a  of  the  mixture  of  centred  Gaussian  laws  depends  on  the  position  x  of  the 
bright  point  of  the  target. 


Case  no.  2. 

By  giving  the  process  with  distribution 


(II. 5)  L  =  l  C(x  )  Sx 

j  j 


we  fix  both  the  laws  of  the  x. 
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and  of  the  C. 


Let  us  take  for  L  a  second-order 
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Levy  process  (cf.  de  Brucq,  part  H) . 

We  introduce  a  positive  bounded  measure  (x,a)  on  M  x  ]R+  satisfying 

2 

the  inequality  /  a  y(dx,da)  <  <=°. 

MxR+ 

By  hypothesis,  for  all  square  integrable  functions  <f> 

2  2 

(i.e.,  /  <j>(x)  a  y(dx,da)  <  °°)  the  second  characteristic  function  of  the 

Mx]R+ 

random  variable 


(II. 6)  <L,<j>>  =  f  (p(x)dL(x)  is  given  by 


(II. 7) 


41 


log  E(exp  i<L,<{>>)  =  / (eaa<^x-* -l-ia<|) (x) )dy  (x,a)  . 


The  stationary  case  corresponds  to  a  product  measure  of  the  type 
dy(x,a)  =  dx  ®  y(da),  and  in  the  Poisson  case  dy(x,a)  =  dx  ®  A6^(da). 

Another  particular  case  is  to  take  the  bright  points  x^  according  to  a 
Poisson  law  of  intensity  A(dx)  on  M  =  [-£,+&]  and  the  amplitude  C  according  to 
a  spherically  invariant  law  with  associated  measure  y(x,dc)  :  the  L6vy  process 
is  here  restricted  to  a  spherically  invariant  mixture.  Thus 
L  =  l  C(x_.)6x_.  must  be  a  Levy  process  with 

^  =  /(eia^x)-l-ia<}>(x)dy(x,a))  = /(e_l-C  ^ 2 -  1) A (dx)y (x,da) ) 

where  the  last  expression  results  from  a  direct  evaluation  of  the  second 
characteristic  function  (of  de  Brucq  and  Gualtierotti ,  p.13).  Replacing 
<j>  by  U(j>  (where  udRl  and  observing  that  /  (eaot0-  1  -  ia0)-^  e"02/2de 
-a2/2 

=  e  -  1,  it  follows  that 

-(u2c2<}>(n)2)/2 


/(« 


-  1)  A(dx)y(x,dc) 


Ja (dx)y (x,dc)  /{eiuCl^X^0  -  1  -  iuc<|>(x)e}  e  6  ^d0 


J  jeiu(c6^  (x) 


1  -  iu(c0)<J>(x)  }A(dx)  y=  e  0  ^dOy  (x,dc)  . 


Making  the  change  of  variables  (c,0)  a  4  c0,  and  putting  y(x,da)  for  the 
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1  -  0  2  /  2 

image  of  the  transition  probability  e  d0y(x,dc)  under  the  change  of 
variables,  we  have  for  a  test  function  <j> 

(II. 8)  ^(u)  =  /(e”U  C  -  1)  X  (dx)y  (x,dc) 

=  J|eiua^*'X'*  -  1  -  iua<j>(x)}A(dx)y(x,da)  . 

This  second  characteristic  function  ¥  (u)  is  remarkable  in  that  it  is 
real  and  depends  only  on  u2.  In  this  case  T  is  invariant  when  u  is  replaced 
by  -u;  so  also  is  the  characteristic  function: 

E(eaU<^’^>)  =  exp  T(u)  =  exp  YC-u)  =  E(e  aU  ). 

The  moments  of  odd  order,  where  they  exist,  are  therefore  zero.  If  this  last 
fact  is  the  only  hypothesis  imposed  on  the  process  L,  we  have  for  all  test 
functions  <J>  and  for  all  ue(E 

Y  (u)  =  J|eiua^  *-X-*  -  1  -  iua<)>  (x)  }dy  (a,x) 

=  J'{e-iua4)  _  i  +  iuat()  (x)  }dy  (a,x)  =  ^(-u) 

=  J { elu("a) ^  ^  _  i  _  iu(-a)cj)  (x)  }dy  (a,x)  . 

In  this  case  the  probability  y  and  the  image  of  y  when  a  is  replaced  by 
-a  lead  to  the  same  second  characteristic  function.  Taking  the  average  of 
these  two  probabilities  we  can  therefore  suppose  that  y  is  invariant  for  the 
change  of  variables  a  -*■  -a  :  the  probability  y  is  symmetric  in  a.  In  this 
case  we  talk  of  a  symmetric  Ldvy  process. 

Case  no.  3 

We  introduce  here  a  generalization  of  Gaussian  measures  which  we  can 
talk  of  as  Spherically  Invariant  Measure. 

Supposing  that  J |G(t ,x) | 2a2v (dx)y (da)  <  °°  we  give  a  sense  to  the  expression 

(II. 9)  Z (t)  =  /  G(t,x)C(dx)  . 

[-4, +4] 
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Thus,  the  points  x  of  the  target  have  a  continuous  instead  of  a  discrete 
random  effect. 

We  introduce,  for  every  function  <fi  of  L2(M,M,v),  a  random  variable  CO), 
centred  with  variance 

(II.  10)  E (C O) 2)  =  /<j>2v(dx)  /  a2y(da) 

M  R+ 

and  with  characteristic  function 

(II.  11)  E (exp  iC O) )  =  /  exp  (—  v y  (da) 

]R+ 

In  the  application  M  =  [-£,+£]  and  M  =  [-& ,  +&] ,  so  that 

(11.12)  Z(t)  =  C[G(t,x) ]  =  /  G(t,x)C(dx)  , 

[-*,+*] 

the  process  4>  e  L2  ([-£,+£],  [-&,+&] ,  v  (dx) )  ->•  C(<j))  is  the  spherically  invariant 
measure;  the  various  points  T  of  the  target  reemit  in  a  correlated  way 
according  to  a  Gaussian  law  with  random  variance  a2. 

The  various  receivers  at  point  S  give  an  observation  Z(t)  which  may  be 
a  scalar  or  a  vector.  In  general  the  random  numerical  values  are  time-averages. 
We  may  write 

(11.13)  <Z,f>  =  /z(t)f(t)dt  where  f  is  a  weighting  function.  The 

observation  Z  thus  becomes  a  linear  process  f  ->  <Z,f>  =  £  C(x.)/G(t,x. )f (t)dt; 

i  3 

the  parameter  f  may  be  for  example  a  square-integrable  function 
(J|f(t)|2dt  <  °°).  The  Hilbert  space  of  the  parameters  f  will  be  called 
H,  in  what  follows. 
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HI .  Variance  and  characteristic  function  of  the  process  Z 

Formally  the  various  cases  correspond  to  the  expression 

(111.1)  <Z,f>  =  l  C(x)jG(x. ,t)  f (t)dt . 

j  J  ^ 

We  will  establish  the  variance  (III. 2)  and  the  characteristic  function 
(III. 3),  taking  the  probabilistic  model  of  paragraph  II. 

(111. 2)  a2 (f )  =  E(<Z,f>2) 

(111. 3)  $(f)  =  E (exp  i<Z,f>) 

Due  to  the  linearity  of  the  process  Z,  we  have  for  all  ueR:  uf  e  H  and 
E(exp  iu  <Z,f>  =  E(exp  i<Z,uf>) 

We  suppose  the  function  x  <G,f>  =  <G(t,x) ,f (t)>  to  be  measurable 
in  x,  and  square-integrable  with  respect  to  the  measures  associated  with  x. 
Case  no.  1. 


a2 (f )  =  E  (/C (X)  f(t)dt)2  =  Jv(dx)  E [C (X) 2 | X-x] <G (t ,x) , f (t) >2 

m  - 

(III. 4)  =  /v (dx)  y (x,da)  a2<G,f>2  . 


Jv(dx)  E  exp  i  C(X)  /  G(t,X)  f(t)dt|X=x 


-a2<G,f>2/2 


=  Jv (dx)  y (x,da)  exp 


Putting  y (dx,da)  =  v(dx)y(x,da)  and  writing  uf  for  f 

(III. 5)  $(f)  =  /  y (dx,da)  exp  :: - j - 2 -  . 

We  have  used  the  hypothesis  of  a  random  variable  C(X)  spherically  invariant 
with  probability  y(x,da)  when  X=x. 

Case  no.  2. 

By  hypothesis,  the  observation  Z  for  the  test  function  f  is  written 


<Z,f>  =  /<G(t,x)  ,f(t)>dL(x)  and  we  suppose  that  /<G (t ,x)  ,f (t) >2  a2  dy(x,a)  < 
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It  follows  that 


(111. 6)  o2(f)  =  J<G(t,x) ,f (t)>2  a2  dp(x,a),  from  the  second  order  expansion 
of  the  second  characteristic  function 

(111. 7)  y(<Z,f>)  =  /(eia<G,f>-  1  -  ia<G ,  f>)dp  (x,a)  . 

This  expression  (III. 7)  is  not  essentially  changed  in  the  special  case 
where  p  is  symmetric  with  respect  to  the  variable  a.  Nevertheless,  it  is 
possible  to  carry  out  the  integration  in  a  on  the  interval  LO>+0°]  instead  of 

On  [-oo5+oo]  . 

A  more  important  special  case  is  the  spherically  invariant  case.  According 
to  (II. 8),  the  second  characteristic  function  is  then 


T  (<Z,f>)  =  /  eia<G»f>_  1  -  ia<G,f>  dp(x,a) 

,  _r2<n  f>2 

=  J(e  ’  -  1)  X(dx)  y(x,dc) , 

from  which  we  obtain,  replacing  f  by  uf 

»<z  f>(u)  =  J(e-W2c2<G,f>2)/2_  1}  i(dx)Y(x,(ic). 

Case  no.  3. 

The  calculation  is  now  carried  out  putting  (x)  =  <G(t,x) ,f (t)>  in 
equations  (I I. 10)  and  (I I. 11).  We  suppose  that  the  function 
x  <G (t , x) , f (t) >  is  square-integrable  on  [-£,+£].  Hence 
/<G(t,x) ,f (t)>  C (dx)  =  C [<G (t , x) , f (t) >]  is  a  centred  spherically  invariant 
random  variable  with  associated  probability  p(da): 


(III. 10)  a2 (f)  4  E 


(J<G,f>  C(dx))2 


J  a2<G,f>2  v(dx)p(da) 


$ (f)  =  E  exp  i  J<G,f>  C(dx)|  =  /  exp  -a2/2  J<G,f>2  v(dx)  pfda) 
and  $£.(u)  =  E[exp  iu<Z,f>3  =  J  exp  F-  U^a  /< G,f>2  v(dx)J  p(da). 
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IV .  Identification. 

It  is  remarkable  that  the  three  cases  considered  lead  to  the  same 
covariance.  For  all  f,g  e  H  we  have 

(IV. 1)  r  (f  ,g)  =  E (<Z ,  f><Z, g>)  =  /  y  (dx,da)  a2<G  (t  ,x)  ,  f  (t) ><G  (t  ,x)  ,g (t)  > 
(for  the  case  no.  3,  y(dx3da)  =  v(dx)y(da).  If  for  all  x  e  [-£,+£]  there 
exists  a  vector  feH  such  that  <G(t ,x)  ,f  (t)>  4  /G(t,x)  f(t)dt  =  0,  then 
<Z3f>  =  0  almost  surely,  for  E(<Z,f>2)  =  0.  Thus  for 

-2y  (R  +  x  sin  0 ) 

G  (t ,  x)  =  P  - - —  cos 

(R  +  x  sin  0 ) 

every  function  f  which  satisfies 

J  cos(ut)  f(t)dt  =  /  sin(mt)  f(t)dt  =  0 
also  satisfies  the  condition  and  so  <Z,f>  =0. 

Hypotheses  of  a  physical  origin  on  the  function  G(t,x)  may  thus  be 
tested  in  this  way. 

The  physical  analysis  of  the  phenomena  establishes  the  form  of  the 
function  G(t,x).  This  function  does  not  appear  to  me  to  be  identifiable 
from  the  covariance  T,  unless  supplementary  hypotheses  are  added. 

The  problem  is  to  identify  the  most  precise  model  compatible  with  the 
experimental  results. 

To  fix  our  ideas,  let  us  start  with  Case  no.  2,  for  which  the  odd 
moments  are  not  a  priori  zero. 

Recall  the  second  characteristic  function 
(III. 7)  T(f)  =  /  [exp  ia<G,f>  -  1  -  ia<G, f>]dy (x,a) . 

Considering,  for  an  ueR,  the  function  uf,  this  leads  to 

(IV. 2)  T^fu)  =  log  E(exp  iu<Z,f>)  =  /[exp  iua<G,f>  -  1  -  iua<G , f>]dy (a,x) . 
The  function  u  ■+  f (u)  can  thus  be  established  statistically  for  a  certain 


2mx  sin  0  2coR  . 
co  t  - - —  + 
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number  of  points  .  .  .  ,  u^.  Making  the  change  of  variables 

(IV. 3)  b  =  a<G(t,x) ,f (t)>,  and  writing  o£  for  the  measure  image  of  y, 

/oo 

o^(db)[exp  iub  -  1  -  iub] 

X  —  oo  t 

For  b  =  0,  the  expression  exp  iub-l-iub  is  zero,  and  so  ck({0})  is  undefined. 

From  a  theoretical  point  of  view,  for  u  e  ]0,°°[  we  have 

+V1  (u)  =  -fC°bzaf (db) exp  iub. 
r  ~oo  r 

The  measure  b2af(db)  is  determined  by  its  Fourier  transform  ^(u)  which 
we  suppose  known.  Thus  the  measure  is  defined  except  for  {0}. 

From  a  practical  point  of  view,  a£(db)  may  be  determined  by  the  method 
of  moments  from  a  series  expansion  in  terms  of  u  (of  Akheizer,  p.22).  Let 
us  calculate  the  moments: 

00  r  •  ~vn 

(IV.  5)  Vf(u)  =  l  -kijL-  Tbna£(db) 


The  estimation  of  the  various  moments  of  <Z,f>  permits  a  limited  series 
expansion  of  V£(u)  in  terms  of  u.  Since  V(u)  =  log  $is  the  second 
characteristic  function  of  <Z,f>,  the  equation 
(IV. 6)  V' (u)  $(u)  =  $' (u) 

gives  us  a  recurrence  relation  for  calculating  the  terms  of  the  expansion 


of  V.  Let 


1U 


.. 


(IV.  7)  *(u)  g  1  +  jt  m.  ■  2!  ‘“2 


t  (iu)z  (iu) 

!  +  — —  nu  +  .  .  .  +  +  *  *  * 


(IV. 8)  ?(u)  =  IT  K1  + 


liHl  k  + 

21  K2 


+  K  + 

+  ~TT~  Ki  + 


In  this  case  the  derivatives  $  *  and  ¥!  are  easily  calculated  and  the 
equation  ¥!(u)$(u)  =  $ 1 (u)  leads  to 


oo  .  p  p- 1 


I 


1*  u1 


K 


=  1  ?  at 


r  i%q 
l  "77 T-fflr 


00  .  £  i- 1 

v  l  u 

=  l  tx~t  mn 


q=0 


«,=! 


(*-1)1 
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or  indeed  to  the  equation  for  the  coefficients  of  powers  of  u 

m„ 


K  mq 

l 

p+q=£ 


(p-l)!q!  («--!)! 


from  which 


V'  £-1 

=  I  C  ,  K  Hi, 

l  L  p-1  p  £  -p 
p=l  r  r  r 


(c11  is  a  binomial  coefficient) 
r 


We  deduce  the  recurrence  relation  for  the  ,  £=1,2,.... 


£-1 


(IV.9)  x  -  m t  -  I  . 

P=1 

Thus  the  moments  m  ,m2 , . . . ,m£  of  <z,f>  allow  us  to  calculate  the 
coefficients  K  ,K, ... ,K  of  ¥ (u) .  The  series  expansion  of  the  two  sides 

Z.  X/ 

of  equation 

(IV. 4)  Vf(u)  =  /  (exp  iub  -  1  -  iub)  a  (db ) 
leads  to  the  equations 
K0  -  0 


K  m. 


Joo 

bo^(db)  by  hypothesis 

i.  -  OO  t 

•  OO  9 

K  =  1  =  j  b  a^(db)  by  normalization 

L  _  oo  t 


Kz  =  =  f°  b*o£(db) 


Proposition  IV. 1:  Assuming  z  arises  from  a  Ldvy  process,  the  matrix 
m(i,j)  =  ^i+j  ’  i > j=l >2, . . . ,n  must  necessarily  be  of  positive  type. 

Proof:  In  fact  for  all  complex  numbers  c^ ,c2> . . ■ ,  c^,  we  must  have 

—  n  i  2 

7  c.c.m(i,j)  =  /  y  c.b1  ap(b)  >  0  .  The  series  expansion  may  seem 
i,j  1 3  °°  i=i  1  1 

a  complicated  detour:  in  fact  it  is  possible,  starting  from  the  equation 
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,  .  ,  ...  iu<z,f>  r00  ,exp  iub  -  1  -  iub.,,  2_  ,,,  , 

(IV. 4)  V  =  log  E(e  )  =  J  ( — - 5 - )t>  a  (db) 

b  • 

to  make  direct  estimates  of  E(e^"U^  Z’^  ^  and  hence  also  of  V^.(Uj)  and  its 
derivatives  V£(up  ,  j=l,2,...,r.  Since  the  observation  <z,f>  is  by  hypothesis 
normalized  (i.e.,  E(<z,f>)  =  0,  E(<z,f>2)  =  1)  we  obtain  K^=0  from  (IV. 9) 
and  we  may  seek  a  probability  b^o  on 3R\{0}  satisfying  /bnb^o^(db)  =  K(n+2). 

-oo  2 

With  these  supplementary  hypotheses,  /  exp(iub)  b  cr~(db)  =  -V'lfu)  and 

—  oo  T  J- 

2 

the  characteristic  function  of  b  is  known  for  values  u^ ,u^ , . . . ,u^ . 

We  have  proved 

Proposition  IV. 2:  Assuming  z  is  due  to  a  Ldvy  process,  and  for  all  f, 
the  function  -V£(u) ,  ueR  must  necessarily  be  a  characteristic  function. 

2 

In  this  case  b  a^(db)  plays  the  role  of  a  spectral  measure  to  be 
identified  from  the  "complex  covariance"  ->i"I(u)  !  If  we  suppose  also  that 
is  absolutely  continuous  with  respect  to  Lebesque  measure,  then 

2  daf  1  r°° 

b  ^  (b)  =  27  L  6XP  iUb  l"Vf(u)]  dU' 

We  have  to  calculate  the  Fourier  inverse  transform  of  the  function 

-^"(u)  known  at  u^ ,u^ , . . . ,u^ .  Rapid  Fourier  transform  methods  may  be  applied. 

We  return  to  the  general  case  of  a  normalized  vector  observation  <z,f  >, 

<z,f  >,..., <z,f  >  associated  with  the  n  test  functions  f  ,f0, . . . ,f  .  The 
’  2  n  1  2  n 

second  characteristic  function,  for  u,,^, . . .  ,u  ‘ 

Y(ulfu2,  ..  .,un)  -=  log  E[exp  i(u1<G,£1>+. .  .+un<G,fn>)  ] 

=  /{exp  ia(u^<G,£^>+ . . •+un<G,fn>)  -  1  -  ia(u^<G,f^>+ . . .+u^<G,f^>}y (dx,da) 

=  /{exp  ia(u10  +...+U  6  )  -  1  -  ia(u. + . . . +u  6  ) }v(d0 ,da) 

Jr  ll  nn  11  nn 

after  the  change  of  variable  x  +  0^  =  <G(t,x) ,f ^ (t)>,  j=l,2...,n  and 
writing  v  for  the  probability  image  of  y. 
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The  probability  v  is  singular  relative  to  the  variables  0^,  j=l,2, . . . ,n, 
and  is  situated  on  the  curve  B:x  0j  (x) ;  thus  only  the  cone  C(0,B)  with 
vertex  0  and  section  B  carries  the  probability 


The  cone  C(0,B)  may  not  be  determined 
from  the  function  T (u^ yu^ , . . . ,u^) ; 
on  the  contrary,  the  deterministic 
physics  of  the  phenomena  must  be  used 
to  define  the  functions  f ^ ,  j=l,2,...,n, 
and  the  response  G.  Hence  the  functions 
fx | — >  G(t,x),fj(t)  are  evaluated  and 
we  deduce  the  curve  B  and  the  cone  (0,B) 
on  which  the  probability  v  must  lie. 


Let  us  further  restrict  the  measures,  in  order  to  arrive  at  the  spheric¬ 
ally  invariant  laws. 

Firstly,  it  is  necessary  that  the  distributions  be  symmetric  about  0, 
i.e.  the  odd  moments  are  zero.  The  equation  T(u)  =  log$(u)  shows  that  for 
replacement  of  u  by  -u,  the  invariance  of  $  is  equivalent  to  that  of  V.  By 

co  2n 

expanding  the  formula  (III. 8)  'F<g  ^(u)  =  I  (-if  7TTT  2/a^n<Gjf>^ndii (x,a) 

’  n=l  ^  '  ' 

2n  2n 

and  by  comparing  with  (IV. 8)  the  moments  satisfy  J a  <G,f>  dy(x,a)  = 

and  proposition  IV  does  not  lead  to  a  new  condition  on  the  model  whereas  in 

the  spherically  invariant  case  new  conditions  appear.  Using  now  formula  (III. 9) 
oo  2n 

T  f  (u)  =  £  (~l)n  —  r  /c^n<G,f>^nX(dx)y(x,dc)  and  comparing  with  (IV. 8) 

<Z,t>  n=l  2n  n* 

we  have 


(IV. 10)  Jc2n<Z,f>2nA(dx)Y(x,dc)  =  K 


K 


2n 


(2n)  !  2n  (2n-l) (2n-3) . .  .3*1 
Proposition  IV. 3  For  a  symmetric  Levy  process  to  be  restricted  to  the  spher- 


K 


2n 


ically  invariant  case,  the  coefficients  — rr— - — r  (instead  of  K„  ) 

(2n-l) (2n-3) . . .3*1  v  2n' 


must  be  the  even  moments  of  a  measure  for  all  neN. 
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Restricting  again  by  supposing  that  the  measure  y  is  of  product  type, 
i.e.  in  the  spherically  invariant  case  y(x,dc)  =  y(dc).  In  this  case 

Jc2n<G,f>2nX(dx)y(x,dc)  =  / c2ny(dc)  /<G , f >2nX (dx)  =  K2n  • 

The  two  measures  y  and  X  are  to  be  identified,  and  this  problem  is  clearly 
undetermined . 

If  we  suppose  that  the  function  x  | - *■  <a(tcx) ,f (t)>  is  known  and 

limiting  ourselves  to  a  stationary  Poisson  law,  intensity  Xdx,  it  follows  that 

(IV. 11)  X  /c2nY(dc)  =  2  y  IV  - 1— —  ;  this  case  is  treated  in 

'  On)!  2n  j<G  f>2ndx 

de  Brueg  §  Gualtierotti . 

We  return  to  the  three  cases  considered  earlier  (sect.  III).  In  (III. 5), 
by  changing  variables  u| - *■  v=u  ;  (a,x)  | - >  a<G,f>  it  follows  that 


(IV. 12)  E(exp  in<Z,f>)=  /exp (-vb) v  (db) ,  where  v  is  the  image  of  y. 

Similarly  in  (III. 9)  we  put  V  =  u2  and  b  =  c2<G,f>  . 

(IV. 13)  log  E(exp  in<Z,f>)  =  / [exp (-vb) -1]  v(db)  with  v  =  measure  image 

of  X*y. 


Also  in  (III. 5)  putting  v  =  u  and  b  =  /<G,f>zv (dx) , 

/CO 

exp(-vb)  v(db). 

Thus  in  all  cases  the  problem  is  to  identify  v  from  the  Laplace  transform 
known  for  a  finite  number  of  points  V2,-#,,Vn  * 

By  a  judicious  choice  of  points  u^,  u2,"’’,Un  We  are  returnec*  t0 
problem  of  moments: 

2  .  .2 

Take  u^  arbitrary;  and  u^  according  to  uj  =  vj  =  •  Thus 


A  a 


E(exp  in  j<Z,f>  =  /[exp  -v^b]^  v(db).  The  change  of  variables 

b  | - >  c  =  expf-v-^b)  leads  to  the  moment  problem  for  the  probability  image 

a:  E(exp  in  j<Z,f>  =  /c^afdc). 
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